Perfect dc conductance of a finite width Mott insulator sandwiched between metallic 
leads at zero temperature: a quantum emergent phenomenon in strongly correlated 

multilayers 
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Using inhomogeneous dynamical mean-field theory, we argue that the normal-metal proximity 
effect forces any finite number of "barrier" planes that are described by the (paramagnetic) Hubbard 
model and sandwiched between semi-infinite metallic leads to always be a Fermi liquid at T — 0. 
This then implies that the inhomogeneous system restores lattice periodicity at zero frequency, has 
a well-defined Fermi surface, and should display perfect (ballistic) conductivity or "transparency" . 
These results are, however, fragile with respect to finite uj, V , T, disorder, or magnetism, all of which 
restore the expected quantum tunneling regime through a finite-width Mott insulator. Our formal 
results are complemented by numerical renormalization group studies on small thickness barriers 
that illustrate under what circumstances this behavior might be seen in real experimental systems. 



In recent years, there has been significant activity in 
the study of multilayered heterostructures constructed 
out of strongly correlated materials. One of the reasons 
for interest in these systems is that new physical phe- 
nomena emerge in the heterostructures that are absent 
in the bulk materials they are composed from, like the 
appearance of two-dimensional electron gases at the in- 
terfaces between band and Mott insulators [1] and their 
low-temperature superconductivity [2] . Here we show the 
possibility of another such striking, quantum emergent 
phenomenon, that at sufficiently low temperatures, a fi- 
nite width Mott insulator sandwiched between two metal- 
lic leads is always described by a robust Fermi liquid fixed 
point, and displays perfect dc conductance. 

While this result may sound counterintuitive, it arises 
in a natural way from the normal-state proximity ef- 
fect. Consider a multilayered device composed of metal- 
lic leads sandwiching a finite number of Mott-insulating 
'barrier' planes. At sufficiently high temperatures, the 
Mott insulator develops a thermal excitation induced 
(one electron) density of states (DOS) within the Mott 
gap, and the system conducts electricity. As the tem- 
perature is lowered, this DOS decreases, but can never 
become smaller than the DOS induced via tunneling due 
to the normal-state proximity effect from the metallic 
leads. The tunneling DOS at zero frequency is expected 
to decrease exponentially with thickness as one moves 
deeper into the barrier planes. But in no circumstance 
does it vanish. Hence, as the temperature is lowered, 
in a pure system, a low temperature Fermi liquid must 
form. Eventually, at T = 0, the electronic self-energy 
vanishes at u) = 0, and the system restores full periodic- 
ity, with a well defined Fermi surface, and no scattering. 
The device develops a perfect conducting channel. The 
temperature scale that governs this behavior is the effec- 



tive Fermi temperature of the innermost barrier plane; 
hence this state is very fragile with respect to increas- 
ing T, to the passage of finite frequency currents, and 
to the presence of a voltage across the planes, disorder, 
or magnetic order in the system (we assume a param- 
agnetic phase in this work). However, for thin enough 
barriers, where the low-energy scale is not exponentially 
small, one can create unique devices with highly nonlin- 
ear properties as a consequence. For example, the sys- 
tem will pass dc current up to some maximal current 
with ease, but then rapidly make a transition to a tun- 
neling response as the current gets too large (a current 
limiter); and similarly if the driving fields change with 
too rapid a frequency, or if the system is in the presence 
of a high enough temperature or applied voltage. While 
this result can become one of merely academic interest, 
not relevant to experimental systems, if the temperature 
scales arc too low, our numerical modeling shows that 
there is a wide parameter range where one should be able 
to see this phenomena. Note that the conclusion of our 
work, that a Mott insulator will always develop a ro- 
bust low-temperature conducting channel, is contrary to 
the conclusion reached in recent work on normal-metal- 
Mott-insulating interfaces [3], where the authors state: 
"our results imply that the Mott insulator is de facto 
impenetrable to the metal" . 

We take the barrier of our metal-barrier-metal sand- 
wich to be described by the single-band Hubbard model, 
and the metallic leads by the noninteracting tight- 
binding model. The Hamiltonian for the system is 
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Here the label a indexes the planes, and the label i in- 
dexes sites of the two-dimensional square lattice in each 
plane. The operator c\ aa (c iaa ) creates (destroys) an 
electron of spin a at site i on the plane a. We set the in- 
plane hopping t" to be nearest neighbor only, and equal to 
t, the hopping between planes, so that the lattice struc- 
ture is that of a simple cubic lattice. The interaction 
term is written in a particle-hole symmetric fashion, with 
U a = U for the barrier planes, and zero for the metallic 
planes. Hence the chemical potential /i is equal to zero 
for half filling, both for the half-filled ballistic metal and 
the half-filled Hubbard system. The system, assumed to 
be at half filling, is then particle- hole symmetric; so we 
drop /z hereafter, and also set t = 1. Furthermore, we 
assume a paramagnetic regime, and drop all spin indices 
in the formulas below. 

The many-body problem is solved within a self- 
consistent inhomogeneous dynamical mean-field theory 
(IDMFT) framework, using the method of Potthoff and 



Nolting [4, 5, 6]. That is, we assume that the self-energy 
is local and the same for all the sites in a barrier plane, 
but can vary from plane to plane, i.e., depends on the 
index a. We label the barrier planes with a values from 
1 to N. Thus a > (N + 1) or a < correspond to the 
metallic layers. The local Green's function at plane a is 
found by employing the quantum zipper algorithm, which 
expresses the Green's function in terms of two continued 
fractions, one from the left, L a (e,uj) and one from the 
right, R a (e, uj) as 

G a (uj) = / dep 2d (e)- r - ( VT~5~T~ \ tTT] > 

(2) 

where P2d{t) is the noninteracting DOS for the square- 
lattice and Z a {uj) = u — 2^ a (uj) , with H a (uj), the 
local self-energy for plane a, being non-zero only for 
a = 1,...,7V, i.e., the barrier planes. The left and right 
functions are given by 



L a (e,uj) = Z a (uj) - e - 
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Furthermore, L a (e,uj) = g ld (w — e) = R a i(e,u>) for 
a < (L a ) and a' > (N + 1) (R a >), where g^M = 
\w ± \^u 2 — 4 is the inverse local Green's function for 
a 1-d tight-binding band. The ± sign in the above equa- 
tions has to be chosen so that in the upper half of the 
complex uj plane L a , R a and g^J are analytic, have a 
positive imaginary part and go asymptotically as u) for 
large \u)\. Thus the IDMFT formalism leads to N sepa- 
rate impurity problems, one for each barrier plane, which 
need to be solved to obtain S Q (w). These are coupled via 
Eqs. (2-4) to determine G a (uj), and hence the N bath 
Green's functions. The details of the IDMFT algorithm 
are well described elsewhere[5, 6], and we do not repeat 
them here. Note that in the present context G a (oj) for 
the metallic planes do not enter the IDMFT recursions, 
and can be calculated at the end for any a. Within the 



IDMFT framework the above procedure therefore fully 
solves the problem of a finite number of barrier layers 
sandwiched between semi-infinite metallic leads. 

On the basis of these equations we argue below that 
the DOS for each plane in the barrier is bounded from 
below by the tunneling result, hence the normalized hy- 
bridization parameter r oQ for the impurity problem at 
plane a is nonzero for all a, and cannot iterate to zero 
via the self-consistent IDMFT algorithm. This then im- 
plies that there is an effective Fermi temperature Tp a oc 
cxp[— 8{7/7rr 0a ] for each of the Mott-insulating planes, 
and when T is lowered below Tp a , that plane will be in 
a Fermi-liquid state, leading to a vanishing self-energy 
at uj = and T = 0. For a thick barrier, with a large 
U, these temperature scales are exponentials of exponen- 
tially large negative numbers, and are hence exceedingly 
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small, and would never be observed in the laboratory. 
But for thin barriers, with U close to the critical U for 
the Mott transition, the effective Fermi temperatures are 
much higher, and could become accessible experimen- 
tally. 

The tunneling DOS, being determined by the proxim- 
ity effect, gets smaller, and hence the self-energy larger, 
as we move deeper into the barrier. We therefore assume 
that 

\Z n/2 (ukO)\ > |S JV/2±1 (o;«0)|>|E JV/2± 2(w«0)| 
> •■• > |£ w /a±i\r/a(w«0)|. (5) 

Hence we can write 

L a (e, lo) ps Z a (uj) - e + 7^ (e,w) (6) 

and 



Ra(e,u) w Z a (u>) -e + 7^(e,w). 



(7) 



Here j%(e,ui) = l/S a _x(o;) for a > 1, and 7f(e, w) = 
-5id(w - e) for a = 1. Similarly, 7„(e,o;) = 
for a < iV, and 7^(e,w) = — <?i<i(w — e) for a = iV. Sub- 
stituting into Eq. (2), using the result |S a (w ~ 0)| ^> 
\j^' R (e,u>)\, expanding the asymptotic form of the frac- 
tion and performing the integration over e then yields 
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s 0). Performing the 
s 0) for a > 1, 7f = 
7^ = 

0.525z. Using this result, we immediately find that 



with 7 ^ = J de P2d (e)jt R (t^ 
integrals yields 7^ = l/E a _i(o; 
0.525i, 7^ = l/S Q+ i(w « 0) for a < N, and - 



T 0Q = ImG - Q V * 0) = Im( 7 £ + 7 f ), 
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with Goq being the effective medium for the impurity, 
which satisfies Gq^ = G^ 1 +S Q . In general, this result is 
a lower bound for Toq, as the DOS may never achieve this 
value as T is lowered from high temperature; this implies 
that our estimates for Tp a will also be lower bounds. 

The bulk insulator phase of the Hubbard model at half 
filling is characterized by a pole in the self-energy that 
takes the form 
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where r(U), the residue of the pole, vanishes when the 
insulating state is no longer stable, i.e., when U < U c \ ps 
11.4 for the simple cubic lattice solved with DMFT (the 
metallic phase disappears for U > U C 2 ps 13.1). For 
the heterostructure, it is interesting to explore the con- 
sequences of the ansatz 
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because this is the form that arises in an exact solution 
of the impurity problem for the Falicov-Kimball model in 
the insulating phase, and in the Hubbard III approxima- 
tion for the Hubbard model with the To a in Eq. (9). In 
this case, the imaginary part of the denominator never 
vanishes, and instead is forced to be nonzero due to the 
normal-state proximity effect. Using this result and the 
definitions of the 7 coefficients, we obtain, for odd N, 



IWwpsO) > 



a-l 



1 < a < 



'0.525/[r([/)] 
1.05/[r([/)] J - 

0.525/[r(U)] N ~ a ffi<a<JV 



N-l 

= ^±i (12) 



Note that this behavior corresponds to an exponential 
decay of the tunneling DOS at lo = 0, with a correlation 
length given by l/lnr(J7). In the Falicov-Kimball model 
on a simple cubic lattice, one has r(U) = (U 2 — 24) /4 [7]; 
for U = 6, just slightly above the critical U for the Mott 
transition at U c = 2^/6, the correlation length is 0.62 [(>] 
and the analytic approach predicts it to be 0.92. For 
the Hubbard model on a simple cubic lattice, a fit ob- 
tained from numerics (see below) gives r(U) ~ (J7 2 — 
C/ c 2 x )/4- 0.13([/-f/ c i)+ 13Ay/U-U cl - 52.31n(J7/C/ c i). 
This suggests that at T = each barrier plane of 
the metal-Mott-insulator-metal structure will become a 
Fermi liquid, leading thereby to a half-filled, simple cubic 
lattice Fermi surface and perfect (ballistic) conductivity. 
(The conductance will be finite, because the ballistic con- 
ductance of the simple cubic lattice is found by counting 
the number of conducting channels [via the Fermi surface 
area] and multiplying by the quantum of conductance - 
this will look like a "contact conductance".) 
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FIG. 1: (Color online.) Local spectral function at the bar- 
rier plane for a single-barrier-plane multilayer, for T = 0.01, 
and a range of U's (increasing from top to bottom near zero 
frequency). 

Next, as a complement to the above arguments, we 
present our numerical work where we solve the impurity 
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problem on each interacting plane with the numerical 
renormalization group (NRG) [8]. We choose the dis- 
cretization parameter A = 2 and we keep 1200 states at 
each NRG iteration, and calculate the local DOS or one 
particle spectral function on different planes, given by 
A a {uj) = —hnG a {uj)/iT. 

Fig. 1 shows our results for the DOS at the barrier 
plane for the single-plane barrier (N = 1) multilayer, for 
T = 0.01 and different values of U. Note how one can 
clearly see the Fermi liquid coherence peak forming for 
all values of U, but the coherence is not complete for the 
larger U values because the effective Fermi temperature 
of the central plane falls below our chosen low temper- 
ature of T = 0.01 for those U values; in addition, note 
how the Fermi peak is narrower for larger U values for 
the same reason. 

Fig. 2 shows our results for a five-plane (N = 5) device 
at three different temperatures for U = 16. We see that 
as the temperature is lowered the spectral functions on all 
the layers tend to the noninteracting simple cubic DOS 
value at u> = 0, but the temperature at which the Fermi 
peak develops is lower for the deeper planes, and the 
width is significantly narrower. We estimate the Fermi 
temperatures for the different planes to be Tpi ~ 0.015, 
T F2 ps 0.0009 and T F3 ps 0.00025. The high energy fea- 
tures in the DOS remain unchanged as T is lowered. 

In so far as transport is concerned, our numerical tech- 
niques are not yet reliable enough to accurately deter- 
mine the conductance of these devices with a Kubo for- 
mula. Nevertheless, we conjecture that decreasing the 
temperature is tantamount to decreasing the length of 
the barrier as more and more insulating layers close to the 
interface become metallic. We speculate that this would 
result in a strong nonlinear dependence of the dc con- 
ductance on temperature, which will be tunneling like, 
but with sharp increases as T is lowered below the effec- 
tive Fermi temperature of each symmetric pair of planes. 
Finally, when T is lowered below the Fermi temperature 
of the central planes, the system will become a perfect 
dc conductor. This is clearly the case in the infinite- 
dimensional limit where vertex corrections vanish. We 
also expect a very narrow Drude peak at low tempera- 
tures, which should be observable using optical or other 
probes of ac conductivity. However, disorder, which is 
always present in any real device, could in practice sup- 
press this low temperature behavior so that the metallic 
regime may be hard to see in real samples, especially in 
those with large Mott-Hubbard gaps, magnetic order, or 
thick insulating barriers. 

In summary, the main result of this work is that a 
system where a Mott insulating barrier is sandwiched be- 
tween metallic leads becomes a ballistic metal at suffi- 
ciently low temperatures no matter how strong the in- 
teraction U , as long as the barrier thickness is finite. 
This striking result is similar to the zero-bias anomaly 
of the quantum dot problem in the Kondo regime [9]. 
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FIG. 2: (Color online.) Local spectral function of the five- 
barrier-plane multilayer for U = 16, and three different tem- 
peratures: (a) T = 0.01, (b) T = 0.001, and (c) T = 0.0001. 
We show the three inequivalent barrier planes (N — 1 is next 
to the interface, N = 3 is at the center). The insets show 
blow-ups of the low frequency region. 



The particle-hole symmetric system (with the same hop- 
ping for all planes and directions) we are considering here 
has the special feature that at zero temperature and zero 
frequency it becomes translational invariant and equiva- 
lent to the noninteracting problem. We also conjecture 
this Fermi-liquid behavior will hold in cold-atom systems 
trapped in a harmonic trap on an optical lattice. The 
central Mott region will be surrounded by normal metal, 
which should create a proximity effect via tunneling, and 
the entire system will become a Fermi liquid as T — > 0, 
so the Mott phase will disappear at zero frequency. It is 
not clear such low temperatures could ever be achieved 
in the laboratory. 
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